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COMPUTER-PLOTTER STUDIES OF VIBRATIONAL EFFECTS IN 
ION-DIPOLE COLLISIONS: "CLASSICAL TUNNELING" 
by John V. Dugan, Jr., and R. Bruce Can right, Jr. 

Lewis Research Center 

SUMMARY 

Computer -plotter studies of ion -dipole collisions are extended to calculate vibra- 
tional effects on capture cross sections and the formation of ion-molecule complexes. 
Two different dipole moment variations are assumed about the equilibrium separation 
of the oscillator, one of which best approximates the quantum time average of the di- 
pole moment. The capture cross section is quite sensitive to the Gaussian dipole mo - 
ment variation p(x), being sharply reduced for a p(x) which is rapidly varying with 
oscillator separation. However, the nature of the ion-molecule complex is relatively 
insensitive to the oscillator model. In certain collisions, the oscillator becomes con- 
strained so it cannot relax to its equilibrium separation; the authors choose to call this 
’’classical tunneling. " The probability of this phenomenon is a function of the oscil- 
lator constants and the ion-dipole interaction. ’’Tunneling" conditions are discussed 
for the two p(x) variations in terms of the ratio of the maximum ion -dipole interaction 
to a "pseudo -oscillator” potential. The latter term depends on the damping constant a 

/_ 2 2s 

(where p = p Q exp' a x ' and p Q is the "static" dipole moment). 

INTRODUCTION 

Capture cross sections and collision complex lifetimes have been calculated numer- 
ically for ion - polar molecule collisions (refs. 1 to 5). The energy transfer in these 
collisions has been studied using computer -plotter techniques (refs. 4 and 5). (The 
polar target molecule is represented by a rotating dipole imbedded in a polarizable 
sphere. The ion is made to reflect at an intermolecular distance R c to simulate the 
hard repulsive cores of the ion and molecule. ) Capture cross sections have been cal- 
culated to set an upper limit to ion-molecule reaction cross sections and to compare 
the predicted ion energy dependence with mass spectrometric results (refs. 6 to 8). 



The permanent dipole introduces multiple -reflection behavior (ref. 3) correspond- 
ing to the formation of long-lived ion-molecule collision complexes. Figure 1 illustrates 
the mechanism by which multiple reflections occur in terms of a plot of energy against 
ion-molecule separation. The rotating dipole alters the potential for the outward tra- 
jectory, outer turning points are introduced, and multiple reflections occur. For the 
Langevin case the ion starts in from infinity (where the pair has relative translational 
energy e), passes over the effective potential barrier at R*, reflects off the barrier, 
and returns to large distances. In a simplified case for the dipole target (i. e. , the 
dipole rotates only slightly), the ion passes over the barrier at R = R f before reflec- 
tion. However, as reflection occurs the dipole can alter its orientation so a new and 
higher barrier appears at a shifted R = R”. The translational energy is now insuffi- 
cient for the ion to escape, and reflection occurs. This discussion neglects energy 
transfer to the rotator, which can lower e and cause reflection even for a fixed dipole 
orientation. 

Ion -molecule collision complexes may be of importance as clustering sites in gas 
discharges or irradiated gases. Such multiple -reflection collisions have calculated 
lifetimes as large as several hundred times a single -reflection period (ref. 4). Motion 
pictures have also been made with the computer to study hindered rotation of the dipole 
by the incident ion (ref. 5). Hindering may favor a specific chemical reaction. 

Lifetimes for multiple -reflection collisions have previously been calculated, in- 
cluding rotational and translational degrees of freedom for the ion-molecule collision 
pair. This report discusses numerical studies of vibrational effects on the capture 
cross section, the probability of multiple reflections, and collision lifetimes. The 
vibrational effects are included by approximating the polar target as a rotating oscillator 
whose dipole moment is a function of internuclear separation along the molecular axis. 


EQUATIONS OF MOTION 


The equations of motion for the interaction of an ion with a vibrating polar rotator 
are derived from the Lagrangian 


L = 


T t + T R + T v 


" V v " V int 


( 1 ) 


where T^ is the relative translational energy of the ion -molecule pair, T^ is the ro- 
tational energy of the polar molecule, T y is the kinetic energy of the oscillator, is 
the potential energy of the oscillator and V int is the ion -molecule interaction term. 
The first term is the relative energy in the center -of -mass (c. m. ) system of ion and 
molecule and can be written 


2 



( 2 ) 


T, = — (X 2 + Y 2 + Z 2 ) 

1 2 

2 2 2 2 

where the ion-molecule separation R = X + Y + Z and m is the reduced mass of 
the ion-molecule pair. The coordinate system is shown in figure 2. The second term is 

T R =-(k 2 + V 2 sin 2 £) (3) 

2 

for a linear molecule, where I is the moment of inertia. The angles £ and rj are 

equivalent to polar and azimuthal angles which define the orientation of the negative end 

of the dipole with respect to a coordinate system fixed on the polar molecule (see fig. 2). 

The kinetic energy of the oscillator T y is defined relative to its equilibrium separation 

r via the reduced coordinate x by 
e J 

rn 0 

T = — - x 2 (4) 

v 2 

The coordinate x is r - r g where r is the instantaneous oscillator separation. Also, 
m Q is the oscillator reduced mass, which can differ from the mass of the polar mole- 
cule. The oscillator potential energy is 

v,^V ! < 5 > 

where k is the force constant of the oscillator. The ion-molecule interaction poten- 
o 2 

tial (refs. 1 to 5) includes the ion-dipole term -pe cos y/R and the ion-induced dipole 
2 4 

(polarizability) -ae /2R terms. The angle y defines the orientation of the vector R 

with respect to a dipole moment vector p which has a direction as shown in figure 2. 

2 2 , 

The dipole moment p equals p Q exp' a x ', where p Q is the "static' 1 dipole moment 
value. It should be noted that the Lagrangian does not include all possible classical 
electrostatic interactions. The interaction between the opposite charge centers of the 
dipole is not considered explicitly. However, the electrostatic interaction along the 
molecular bond is accounted for implicitly by choices of both the force constant k Q and 
the p(x) variation. This is a semiclassical approximation since it says something about 
the quantum mechanical bonding along the molecular axis. The total interaction potential 
is 
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( 6 ) 


V int = 


/ 2 2x 

M Q exp^" a x 'e cos y ae 2 


R 


2R^ 


The total Lagrangian in the center -of -mass system is 


t m , • 2 '2 I ,' t 2 2 . 2 ^ 1 -2 

L = — (X +Y + Z ) + — (£ + v sin £) + — m x 

2 2 2 ° 


2 2 > 


1 , 2 ^o ex P ( " a X )e cos y ae 2 


- - kx + 
2 0 


R 


+ (7) 

2R 4 


where mR cm /2 (the kinetic energy of the center of mass in the laboratory system) is 

omitted since R =0 and 
cm 

cos y = _ ‘ hL = X sin £ sin ?7 - Y sin | cos 77 + Z cos £ 

Rju R 

The equations of motion are obtained from 


d_ /3L\ _ / SL N 
dt \ 3 V \ 3 q i/ 


= 0 


where qj = X, Y, Z, 77 , and x. These six second-order equations give 12 first- 
order equations which are solved for the velocities and coordinates employing the 
variable -step-size Runge-Kutta integration routine described in reference 9. 

The equations of motion are 


X 


Y 


Z 



(R 2 sin | sin 77 - 3XF y ) - 2ae2x 

mR 6 

(9) 

2 

(-R 2 sin ij cos 7 7 - 3YF y ) - ^ ae Y 

( 10 ) 

mR® 


2 

(R 2 cos i - 3ZF ) - 2ae Z 

( 11 ) 

mR 6 
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with 


F y = (X sin £ sin 77 - Y sin £ cos 77 + Z cos £) 


£ = sin £ cos £ + / (-Z sin £, + X cos f; sin 77 - Y cos i- cos 77) 

\m 3 / 


( 12 ) 


” _ / j*e\ (Y sin 77 + X cos 77) - 277I cos | 

\m 2 / sin « 


(13) 


x = . 2a fte cos y + k \_x_ 

V R 2 7 m o 


(14) 


MOTION-PICTURE MODEL 

Computer -made motion pictures of the oscillator collision were made for collisions 
between equal mass ions and molecules. Thus, in this case, each particle is equidistant 
from the center of mass. A sample frame is shown in figure 3. The ion and polar mole- 
cule models consist of a positive sign and a dipole, each surrounded by dotted circles of 
varying radii. The ellipse at the top of the frame indicates the value of x, that is, the 
difference between the instantaneous oscillator separation and the equilibrium oscillator 
separation. The motion pictures were used to detect the tunneling phenomenon which is 
not evident in the time history plots because the x variation is very rapid at small R. 


OSCILLATOR MODEL AND INITIAL CONDITIONS 

At high frequencies (i. e. , on the order of molecular vibrations and fast rotations) 
the dipole moment of a polar molecule is not constant but is a function of its internuclear 
separation. The separation of interest is usually that along a particular interatomic 
bond in the molecule (ref. 10). The variation in dipole moment p(x)is assumed Gaussian 

( a 2 x 2 i 

for these studies (i. e. , ju = M o exp v ’), where p o is the "static” dipole mo- 
ment, a is the damping constant, and x has been defined. The values of u are 

OQ ° 

3. 92 Debye units (1. 3x10 C-m) for CHgCN and 0. 10 Debye units for CO. Such a 
variation is a reasonably representative approximation for diatomic and polyatomic 
molecules of interest (ref. 10). (The maximum x value for the free oscillator, set by 
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values of k Q and initial oscillator energy E Q , was 0. 05 A (0. 005 nm) in all cases re- 
ported. ) To test the sensitivity of the results to the shape of the p(x) variation, the 

damping coefficient a was selected so that p varied gradually with x for oscillator 
o 

set A. This a value is 160. For oscillator set B, p varies more rapidly with x 

o 2 3 

(p = 0. 02 p Q at x = ±0. 05 A with a = 1. 6x10 ). The shapes of the two Gaussian varia- 
tions are shown in figure 4. 

The most realistic model is that which most nearly simulates the time average of a 

quantum oscillator (i. e. , some P e ff( x ))- The M eff ( x ) is the product of p(x) weighted 

by the probability P(x) of finding the oscillator at a given x value. The P(x) is 

| i// 2 (x) | for a quantum oscillator and this function is sharply Gaussian for a ground-state 

oscillator (ref. 11). The classical probability P(x) weights the turning points heavily. 

The quantum and classical probability functions P(x) are shown in figure 5 (ref. 11). 

2 

The classical time average of p(x) with small a is probably a better approximation 

o 

to the realistic p g ^(x) for ground-state molecules than p(x) with large a . In the 
latter case, small p values are weighted heavily. The oscillator spends much time at 
the turning points. These conclusions apply strictly only to a free oscillator. The 
interaction can change the p g ^(x) considerably. 

The period of change in p is t q « 2tt y m o /k Q , where k Q and m Q are the oscil- 
lator force constant and reduced mass. The values of k Q and m Q for the oscillator 
define the maximum values of x and x, for fixed E 0 . The m Q value is roughly the 
same as the mass of the target molecule and is representative for a diatomic oscillator 
like CO. The ratio k Q /m 0 was chosen so as to reduce computer time; thus, the k Q 
is somewhat smaller than the value of a representative polar molecule. The energy 
values E q chosen were 0. 36 eV and 0. 036 eV, corresponding to stiff and soft (nearly 
thermal) oscillators. However, only 50 cases were studied for the stiff oscillators be- 

4 

cause computer time was prohibitive. A k Q value of 4. 8x10 dynes per centimeter 
(0. 48 N/cm) was assigned for the calculations reported here, with an m Q value equal 
to 10 times the mass of the hydrogen atom. These values correspond to a t q of 
~10‘- 13 second. 

Integration for trajectories of all oscillator cases causes extended computer running 

4 

times, a minimum of 5 minutes for v = 5x10 centimeters per second. The oscillation 
period is only 10” 3 of the characteristic rotational period where = ir y/21/E-^, 

Ej^ is the rotational energy, and I is the moment of inertia. Therefore, a smaller 
average time step is required to perform the integration. The oscillator motion some- 
what separates out from the rotational and translational motion; however, it will be 
shown that the nature of the p(x) variation can be important in determining the ion 
turning points. Individual collisions with stiff targets require as much as 10 minutes 
to 1/2 hour computer time, which is clearly prohibitive since hundreds of cases must 
be studied. 
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The collision results reported herein are for the soft oscillators with E Q = 0. 036 eV. 
The total energy of the free oscillator is 

E = — m x 2 + - k x 2 (15) 

o 2 0 2 ° 

where all quantities have been previously defined. The maximum displacement of the 
free oscillator at large distances from the ion is given by x m = V 2 E oA 0 (*• e - , th e 
separation where the oscillator energy is purely potential with no interaction). The 
initial separations x were chosen randomly in the interval |x| < x m . Approximately 
400 collisions have been studied for impact parameters b from 3 to 12 A. These b 
values correspond roughly to impact parameters for probable capture with CO and 
CHgCN targets (refs. 1 to 3). The oscillator energy remains positive throughout inter- 
action but sometimes decreases in value. Such behavior is, of course, not possible for 
a quantum oscillator in the ground state, which must have energy 1/2 hi/ Q . Thus the 
observed ’’cooling” of such oscillators applies only to molecules which are initially in 
excited vibrational states. The classical oscillator approximation is most adequate in 
describing variations in energy of soft oscillators since quantum gaps are relatively 
small. 

The differential equation of motion for the oscillator is 

^25^ + kA — (14) 

! ) m o 

The first term is the ion-dipole interaction term, whereas the second term is the con- 
ventional restoring force term for the oscillator. 


X= _ [ 2a_pe 


R 


RESULTS 

Multiple-Reflection Behavior 

The probability of multiple reflection is the ratio of the number of multiple - 
reflection collisions to the total number of collisions studied for a constant impact 
parameter. This probability for oscillating polar targets f^ is compared with results 
for purely rotating targets in table I. These results are for an initial ion velocity of 

4 

5x10 centimeters per second, 25 cases per point. In addition to choosing the initial 
value of x randomly, the values of initial positions and energies are chosen randomly 
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(ref. 9). The results of set A are for an a value of 160 which makes p = 2p /3 at 

x . Results are presented for both CO and CHoCN. Results for set B are for an a 
m s o _o 3 

value of 1. 6x10 A which makes p. = 0. 02 p 0 at x = ±0. 05 A, and include only CHgCN 

results. The effect of oscillations clearly depends on impact parameter; however, it is 

difficult to detect general trends in the results. The variation p(x) for the oscillator 

is important; this can be seen from the difference between identical random number sets 

for different a values at roughly the same b values. It should be pointed out that 

although the number of cases studied is limited, the variation between different random 

number sets of 25 cases (without the oscillator) is roughly 15 percent. 

The f-^ value is zero at b = 11 A because the capture probability for set B drops 

off abruptly in this range. This dramatic effect on the .apture cross section is demon- 

2 

strated by the plots of capture ratio Cj^ against the square of the impact parameter b 
for CHgCN in figure 6. The capture ratio is the fraction of collisions in which the ion 
and the molecule approach within a specified separation R (usually several A). The 
capture cross section a c for ion-dipole collisions is simply related to Cj^ by 

r 2 

a =n J C R (b)d(b ). There are fewer points for the vibrating targets, only 25 cases 
c 0 K 

per point compared to 50 cases (for the purely rotating targets) because of the prohibi- 
tive computer running times. The results agree at low-impact parameters b = 6 and 
8 A. However, there are marked differences between results for b = 10 A as well as 
for b = 11 A. The differences between values are well outside the rough 10 percent 
variance in cross section between 25- and 50-case sets of trajectories. The oscillator 
studies were not done for b > 11 A since the C R value approaches zero. The effect 
of the oscillator on capture ratio will be much less for CO than for CHgCN since the 
polarizability term determines the capture cross section for CO (refs. 4 and 5). 

The collision lifetimes are relatively insensitive to vibrations and the average time 
for multiple -reflection collisions is 2x10”'^ to 3x10”^ second for CHgCN targets. It 
should be noted that the trajectories were cut off at a maximum of 10 reflections so that 
no 'meaningful maximum reflection time was calculated. Results of previous studies 
indicate that the average number of reflections is 15 to 25 for certain b values (ref. 5). 

Certain features of the energy exchange between the translational, rotational, and 
vibrational degrees of freedom are shown in the plots of figures 7 and 8. The features 
of the variations in relative velocity and rotational energy in figure 7 are not affected 
by the presence of the oscillator. The modulation of the reduced oscillator separation 
x via the ion -dipole interaction is evident in the region R « 7 to 8 A, as shown in fig- 
ure 8(a) for a CHgCN target. The corresponding variation in oscillator energy is shown 
in figure 8(b). The time interval of the modulation envelope is roughly equal to the 
rotational period. 
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Multiple -reflection behavior is shown in figure 9(a) for a representative collision. 
The hindered rotation accompanying multiple reflection is shown in figure 9(b). The 
symmetrical nature of the net spiraling due to multiple reflections (refs. 2, 3, and 5) 
is somewhat distorted by the oscillator. However, the general features of this phenome- 
non are unchanged by the oscillator, as is hindering of the dipole. The hindering of the 
CHgCN target occurs at large ion-molecule separations (~7 to 8 A), as can be seen from 
the decreasing envelope of the ion-dipole orientation y plot in figure 9(b). The charac- 
teristic ion turning points of the ion-molecule complex are large (10 to 20 A for CHgCN) 
and their values are somewhat insensitive to the presence of the oscillator. 


"Classical Tunneling" 


Vibrational motion has been added to the ion -dipole motion pictures which were 
described in references 5 and 12. The modulation of the oscillator which is suggested 
in the time history plots is confirmed in the motion pictures. This modulation is, of 
course, dependent on the value of cos y (i. e. , the ion-dipole orientation angle which 
depends on the relative motion of ion and molecule). 

In viewing the motion pictures it was noted that there is a change in the mode of the 
oscillation in several instances. This change in mode is such that the oscillator cannot 
relax to its equilibrium position. This behavior is somewhat analogous to quantum 
mechanical tunneling through a potential barrier so it is herein called '’classical tunnel- 
ing. ” This change in mode occurs when the oscillator acceleration goes through zero; 
that is, when the restoring force is balanced by the interaction term. It should be noted 
that x can be zero where tunneling does not occur (i. e. , the free oscillator when x = 0). 

It is useful to discuss tunneling (x = 0) by inspecting the condition that the two terms 
on the right side of equation (14) cancel, that is, 


2a cos y ( k _ Q 
9 O 

Pt 


(16) 


or rewriting 


-% e cos y 2a 2 „. rr (-a 2 x 2 ) . , 


R 


(16a) 
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Taking the natural logarithm of both sides in equation (16a) we obtain the condition 


In 



2 2 n 
-ax = 0 


( 17 ) 


where 


V 


m 


-fX Q e cos y 


R 


2 


(18) 


and 


V 


o 



(19) 


2 2 

Since ax is always positive, it is clear that cos y must be less than 0 for tunneling 

to occur. The logarithmic argument is the ratio of the maximum interaction potential 

V to V , the latter can be viewed as a pseudo-oscillator potential. Tunneling is clearly 

quite sensitive to the value of a , and the range of a x^ is quite different for the two 

sets of oscillators. It is useful to inspect the limits on R and rough constraints on y 

2 2 

for tunneling at the most probable a x values. 

2 2 2 2 

Although 0 < a x <1.5 from set A, the mean a x value for the nearly free 

oscillator is approximately equal to 1 since it spends most of its time near the turning 

2 2 

points. The condition (ax) =1.5 puts an upper limit on R for set A tunneling at 

III ctX 

29 A. At large R where V is approximately equal to V the tunneling condition 

in o 22 

suggests cos y values near -1 for the representative value of a x . At small R values 
the -ion tends to orient the dipole attractively, thus not meeting the tunneling require- 
ment that cos y be less than zero. 

2 2 2 2 

The range of a x values for set B is 0 < a x <4; however, the larger values are 
again favored because of the nature of the weighting function P(x). The upper limit to 
R for tunneling in this case is only 15 A. Since relatively large values of V /V_ 

1 v 111 O 

(where V Q « KT'kT) are allowable for R « 15 A tunneling appears more probable for 
set B. 

The time history plots of oscillator coordinate x and velocity x in figure 10 show 

that x stays near zero for a significant fraction of an oscillation period (where x is 

nearly constant). All tunneling is observed when y ~ n/2 for set B oscillators (large 
o 

a ) so that cos y is very small and negative. No tunneling is observed for targets from 
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oscillator set A, but motion pictures were used to detect tunneling in only five cases. 
It should be pointed out that this phenomenon is less probable than multiple reflection 
and seems to depend critically on the nature of the dipole moment variation. 


CONCLUDING REMARKS 

Preliminary studies of vibrational effects on multiple reflections, collision lifetimes, 
and capture cross sections in ion-dipole collisions have been completed. For thermal 
oscillators (E Q = 0. 036 eV) the effect of the vibrator depends on impact parameter as 
well as on target species. The fraction of multiple -reflection cases is a complicated 
function of the oscillator model, particularly the shape of the dipole moment variation 
p(x). The oscillator can have a dramatic effect on the capture ratio Cj^; this corres- 
ponds to a lowering of the capture cross section, which is particularly sensitive to the 
damping constant a. The inclusion of additional degrees of freedom (e.g. , coupled 
oscillators) might markedly increase the collision complex lifetimes. Vibrational 
studies must be extended to a variation in parameters before more general conclusions 
can be drawn about the role of the oscillator. "Classical tunneling" of the oscillator via 
the ion -dipole interaction suggests a mechanism for rearrangement and/or dissociation 
which operates without the force constant changing. Quantitative conclusions regarding 
this and related phenomena require further study, particularly in assessing the validity 
of the semiclassical oscillator model for real polar targets. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 23, 1971, 

129-02. 
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APPENDIX -SYMBOLS 


a 

damping constant in dipole mo- 
ment expression, A - * 

t r 

rotational kinetic energy of 
polar molecule, eV 

b 

impact parameter, A 

T t 

relative translational kinetic 

C R 

capture ratio, dimensionless 


energy, eV 

E 

oscillator energy, eV 

T 

V 

oscillator kinetic energy, eV 

e 

-IQ 

electronic charge, 1.6x10 v C 

V int 

interaction potential, eV 

f 

multiple -reflection fraction 
dimensionless 

V m 

maximum ion -dipole potential, 
eV 

k 

_Q4 

Planck's constant, 6.623x10 
J-sec 

V o 

pseudo-oscillator potential 
(V 0 = k 0 / 2a2 ), eV 

k 

oscillator force constant, N m - ^ 

V v 

oscillator potential energy, eV 

L 

Lagrangian, eV 

X 

translational coordinate of ion, A 

L’ 

translational angular momentum 
2 -1 

(L* = mvb), kg m sec 

X 

reduced oscillator separation 
(x = r - r e ), A 

m 

reduced mass of ion-molecule 

x m 

maximum value of x, A 


pair, kg 

Y 

translational coordinate of ion, A 

m 

0 

oscillator mass, kg 

Z 

translational coordinate of ion, A 

R 

ion-molecule separation, A 

a 

electronic polarizability of 

R c 

ion-molecule reflection distance, A 


polar molecule, A 

R 

center- of- mass coordinate in 

P 

quantum oscillator coefficient, 

cm 

o 

laboratory system, A 


A” 2 

R T ,R M 

separations corresponding to 

y 

ion -dipole orientation angle, deg 


maxima in effective potential 
for ion-dipole interaction, A 

e 

relative ion-molecule transla- 
tional energy, eV 

R* 

separation corresponding to 
maximum in effective potential 
for Langevin collision, A 

V 

M 

azimuthal angle for rotational 
motion, deg 

instantaneous value of dipole 

r 

instantaneous oscillator separa- 

o 

tion, A 


moment, Debye units 
(1DU = 3. 3xl0“ 30 C-m) 

r e 

equilibrium oscillator separa- 
tion, A 

% 

static value of dipole moment, 
Debye units 
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oscillation frequency, Hz t q 

polar angle for rotational motion, u> Q 

deg 

capture cross section, A 


period of free oscillation, sec 

oscillation frequency 
(u> o = 2ttu q ), rad sec - * 
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TABLE I. - MULTIPLE -REFLECTION RESULTS FOR 


HARMONIC OSCILLATOR STUDY, 25 COLLISIONS 
STUDIED PER IMPACT PARAMETER 


r 4-1 

[Ion velocity, 5x10 cm sec ; oscillator force constant 
k Q = 4. 8x10* N m”*; reduced mass of ion- molecule pair 
m Q = 1.64X10' 26 kg.] 


Polar 

Impact 

Fraction of collisions with 

Static 

molecule 

parameter 

multiple reflection, 

dipole 


A(10 _ * nm) 

f R 

moment, 



Without 

With 

4 0 

Debye units 



oscillator 

oscillator 




Set A (a 2 = 160 A 2 ) 

. 


CO 

3. 0 

0.24 

0. 16 

0. 10 


4.0 

.20 

. 16 


10 


5.0 

.24 

. 16 


10 

CH 3 CN 

5.0 

0.24 

0.20 

3. 

92 


6. 0 

. 16 

. 08 



1 

7. 0 

. 36 

.48 




9. 0 

.28 

.36 




11. o 

.20 

.20 




Set B (a 2 = 1. 6xl0 3 A' 2 ) 



CH 3 CN 

6. 0 

0. 16 

0.60 

3. 

92 


8. 0 

.24 

. 12 




10. 0 

.24 

. 32 




11. 0 

.20 

.0 





CO 



R ► CS-51020 


Figure l. - Schematic diagram illustrating mech- 
anism of multiple-reflection phenomena in ion- 
dipole collisions with plots of effective potentials 
for Langevin and permanent dipole collision 
systems against ion-molecule separation. 



Figure 2. - Coordinate system used for study of ion-dipole interaction 
with oscillating polar targets. 
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RATIO OF DIPOLE MOMENT TO STATIC 



Figure 3. - Sample motion-picture frame for collision study with 
oscillating targets, showing projections of ion and polar mole- 
cule and location of oscillator reduced mass at top of frame. 
{Oscillator separation scaled by factor of 32 for clarity of dis- 
play. ) 



Figured - Dipole moment variation for polar target oscillators studied (sets A and B); |x m | 0 = 0.05 A; 
p - p^expl-a 2 ^). 
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(A) QUANTUM OSCILLATOR; (Jjj - (B) CLASSICAL OSCILLATOR; 

y^n exp(-px^); P(x) = ipQ y55j}; p= P(x) = [rt 2 Eg/kQ - x?) 11 ^ \ 

moa) 0 /2fi, WHERE w 0 IS OSCIL- 
LATOR FREQUENCY IN RAD SEC ' 1 
(REF. 11). 

Figure 5. -Comparison of probability densities for quantum and classical 
ground-state oscillators, as functions of reduced oscillator separation. 



Figure 6 . -Comparison of capture ratio variation with impact parameter for oscillating (from set B) 
and pure rotating polar CFI 3 CN targets. Results are for initial ion velocity of 5x10^ cm sec -1 

where capture cross section a c = * /*“ C R (b 2 )d(b 2 ). 

• / 0 
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ROTATIONAL ENERGY SCALED ABSOLUTE VALUE OF VELOCITY 

TO kT p = 0 026 eV SCALED TO INITIAL VALUE OF 

5. 5x104 CM SEC-1 



< 

Of 

< 



20 15 10 0 5 

ION -MOLECULE SEPARATION, R, A 

(A) OSCILLATOR SEPARATION. 



25 20 15 10 5 0 

ION-MOLECULE SEPARATION, R, A 

Figure 7. - Variation of relative velocity and polar molecule rotational 
energy during CH 3 CN - parent ion multiple-reflection capture col- 
lision with oscillating target 



(B) OSCILLATOR ENERGY. 

Figure 8. - Variations in oscillator separation and oscillator energy before 
reflection in CH 3 CN - parent ion capture collision with oscillating target 
from set B (a2 = 1.6xlCP A“2). 




OSCILLATOR SEPARATION, 



X AXIS 

(A) PROJECTIONS OF ION TRACE 
IN CARTESIAN XZ PLANE OF 
COORDINATE SYSTEM WITH 
MOLECULE AT CENTER. 



(B) ION-DIPOLE ORIENTATION ANGLE. 

CS- 56280 


Figure 9. -Evidence for multiple reflection and hindered rotation behavior 
in CH 3 CN ion -dipole collisions with oscillating target. 



TIME, t, SEC 

CS- 56282 


Figure 10. - Variations of oscillator separation and velxity for ion-dipole collision demon- 
strating "classical tunneling. " 
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